Factorization, Power Corrections, and the Pion Form Factor 
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This paper is an investigation of the pion form factor utilizing recently developed effective field 
theory techniques. The primary results reported are: Both the transition and electromagnetic form 
factors are corrected at order A/Q due to time ordered products which account for deviations of the 
pion from being a state composed purely of highly energetic collinear quarks in the lab frame. The 
usual higher twist wave function corrections contribute only at order A 2 /Q 2 , when the quark mass 
vanishes. In the case of the electromagnetic form factor the A/Q power correction is enhanced by 
a power of l/a 3 (Q) relative to the leading order result of Brodsky and Lepage, if the scale \/AQ is 
non-perturbative. This enhanced correction could explain the discrepancy with the data. 



INTRODUCTION 

Making predictions for hadronic observables is ex- 
tremely difficult given the complexity of the theory of 
strong interactions (QCD). In the absence of a solution 
to the theory, we are forced to accept reduced predictive 
power. We can make predictions only after we have ex- 
tracted some crucial related information from the data. 
Moreover, the necessary information is usually in the 
form of non-local matrix elements and not just fixed 
couplings. Even with these lowered expectations, it is 
still highly non-trivial to find observables which we can 
predict from first principles. The primary tool at our 
disposal is factorizationjlj. Observables which are "fac- 
torizable" can be separated into long and short distance 
contributions. Asymptotic freedom allows for a calcula- 
tion of the short distance piece via a perturbative expan- 
sion in the coupling, while the long distance piece is in 
principle calculable, but in practice must be extracted 
from experiment. The predictive power lies in the fact 
that factorization implies universality. That is, the same 
non-perturbative factor can appear in predictions for dis- 
parate processes. 

Proving that a certain observable is factorizable is a 
highly non-trivial process. The pioneering works on the 
subject 0, IH were based upon diagrammatic tech- 
niques and can be quite intricate. Using these techniques, 
factorization was shown to leading in 1/Q, where Q is 
the large energy scale in the process. For example, the 
photon-pion transition form factor which is defined via 
(q =P* -p-y) 
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can be written to leading order (in 1/Q = 1/ y/\ q 2 |) as 
2U 
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while the pion EM form factor defined as 

(tt(p)I ^(0)|7r(p')) = ^(Q 2 )(PM+^) (3) 
may be written as 
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In these expressions cj>„ carries the universal, non- 
perturbative, information about the structure of the pion 
and is defined via 
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where W is a light-like Wilson line needed for manifest 
gauge invariance. C\ and T\ are the perturbatively calcu- 
lable, process dependent, high energy Wilson coefficients. 
At asymptotically large values of Q, the pion wave func- 
tion is dominated by its first moment and approaches the 
form cj)„(x) oc x(l — x)f„. Experimentally, the prediction 
for the transition form factor agrees at the ten percent 
level, while the EM form factor is off by a factor of order 
one. Thus, the relevant questions is, how large are the 
subleading corrections in each case? This paper under- 
takes the task of systematically categorizing the power 
corrections to these processes, including the contribution 
from the so-called "end-point region" . The issue of the 
breakdown of factorization will also be discussed. 



SCET 

The results in this paper are derived utilizing recent de- 
velopments in the so-called Soft-Collincar Effective The- 
ory (SCET) 6] . In this approach, factorization proofs 
simplify because modes with varying virtualities are sep- 
arated at the level of the Lagrangian 7J . Proving factor- 
ization becomes tantamount to determining if the theory 
properly accounts for the IR physics of the process under 
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consideration. Perhaps the true power of the effective 
field theory approach, as applied to exclusive processes, 
is that it enables one to account for power corrections in 
a systematic fashion including the so-called "end-point" 
contributions 1 . 

To build the proper effective field theory one isolates 
the degrees of freedom responsible for the non-analytic 
behavior of the low energy theory. Usually this cor- 
responds to simply integrating out massive fields, but 
there are cases where one wishes to explicitly sep- 
arate certain subset of fluctuations of the fields. This 
separation allows for manifest power counting, which in 
turn makes the classification of power corrections rela- 
tively simple. In SCET the relevant modes in the 1R are: 
collinearpc = (p+ , p~ , p ^ ) = (np c ,np c ,p^) ~ Q(A 2 ,I,A), 
soft p s ~ Q(A, A, A) and usoft p us ~ Q(A 2 , A 2 , A 2 ), where 
n 2 = n 2 — 0, n-h = 2, and A <C 1 is the expansion pa- 
rameter. The relevant modes are fixed by the external 
momenta and use of the Coleman-Norton theorem. In 
0, the authors chose A cx A/Q, as this fixes the trans- 
verse momenta of the external lines to be of order A, 
as it should be physically. This also means that the 
usoft modes have virtuality less than A, that is, their 
wavelength is longer than the confinement radius. These 
"hyper-confining modes" should not contribute to phys- 
ical processes. In all the cases discussed in Q, it was 
shown that these modes do indeed cancel in the matrix 
elements. There is one additional mode which may con- 
tribute in the case when there are at least two collinear di- 
rections. Namely, the so-called "Glauber region" , whose 
momenta scale as (A 2 , A 2 , A). While these modes typi- 
cally are not relevant in exclusive processes, this has yet 
to been shown within the context of SCET. 



The Need for an Intermediate Theory 

New subtleties in the theory arise when one consid- 
ers power corrections to exclusive processes due to what 
is known as the "end-point" region, which is defined as 
the contribution to the process where the hadron is in 
an asymmetric configuration. That is, one or more of 
the constituents carries a parametrically small part of 
the longitudinal momentum 2 . In the context of effective 
field theory language, this configuration corresponds to 
a sub-leading fluctuation. To see this note that to pre- 
serve manifest power counting, we take the pion to be 
an eigenstates of the leading order Hamiltonian of the 
effective theory, which as will be discussed below, only 



contains collinear fields. As such, the overlap of the pion 
with a field which contains a net soft quantum number 3 
will be nil. Thus, in the effective field theory formalism 
these soft regions of the wave function show up via time 
ordered products wherein a collinear field fluctuates to a 
soft field and then back to a collinear field. These time 
ordered products must involve an operator which cou- 
ples a collinear field with a field whose momenta scale 
as (A, A, A). However, momentum conservation forbids 
the coupling of such a field with a collinear field whose 
momentum scales as (Q, A 2 /Q l A). 

This issue was addressed in where it was pointed 
out that such interactions will be properly accounted 
for, if one considers working in two stages. In the ini- 
tial theory, labeled SCETI, the scaling parameter A is 
order yj A/Q, this theory is valid at scales below Q but 
above \/AQ. In this theory interactions between collinear 
(1,A 2 ,A) and ultra-soft (A 2 , A 2 , A 2 ) gluons are permit- 
ted. At the scale y/AQ, we match onto a second the- 
ory SCETII, where A now scales as A/Q. In doing this 
matching we integrate out modes with invariant mass 
q 2 > AQ. In SCETII, there are no interactions between 
collinear and soft modes 4 



THE TRANSITION FORM FACTOR 
Leading Order 

In the effective theory each mode is interpolated by 
a distinct field and scales homogeneously in A. For in- 
stance, collinear modes with large light-cone momentum 
in the n direction are interpolated by £ n and A£ for 
fcrmions and gauge bosons respectively. These fields have 
support only over momenta of order A 2 , as their large 
light cone and transverse momenta have been scaled out. 
The leading order Lagrangian for these collinear fields 
with momenta in the n direction, is given by 
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The operators V and V ± are derivative like operators 
whose eigenvalues are the large light cone and transverse 



1 These contributions are sometimes also referred to as the "Feyn- 
man" or "soft" regions. 

2 I will always be discussing processes/frames where the hadron is 
nearly light-like. 



3 That is, an operator which has no overlap with a purely collinear 
state. 

4 One can think of the ultra-soft modes in SCETI as becoming the 
soft modes in SCETII. 
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momenta respectively. The gluon action, ci°J , can be 
found in Q. 

Let us start by considering the pion transition form 
factor. The leading order matching for this process was 
performed in [tJ. This process involves the scattering 
of a highly virtual photon and a quark-anti-quark con- 
stituent pair of an on-shell photon in the n direction. We 
will work to leading order in a em . 5 By integrating out 
the hard off-shell intermediate state in the 7* + 7 — * q + q 
process, we generate a two quark operator in the SCETI. 
The most general spin structures for currents with two 
collinear particles moving in the same or opposite direc- 
tions are [jj 



£n {$,$75 ,$7xH« 5 Cn{l,75,7x}^ 



(8) 



From this result we can see that, for the case at hand, 
there is only one relevant operator structure which inter- 
polates for the pion namely $75, thus the leading order 
matching result is of the form 



(9) 



where the isospin structure has been suppressed, and 
Here the Ws are the Fourier trans- 



forms of light-like Wilson lines 



W(y, 00) = Pexp ( i / n ■ A n (Xn)dX 



(10) 



and u>i : 2 the total collinear momenta of the jet like struc- 
ture Wt£ and £W respectively. Typically each operator is 
accompanied by label subscripts, such labels are implied 
if they are not explicit. There is an implied sum over all 
label momenta such as ^1,2, which will be restricted by 
momentum conservation when we take a matrix element. 
We may then decouple the usoft modes from the collinear 
modes in the Lagrangian via field redefinitions [5] . 

i(n,n) — ^( n ,R)^(n.«) ' ^{n,n) = ^"( n ,n)^(n,n)^(n.n) >i^) 

where Y is an usoft Wilson line defined as 



Y„ = Pexp ( i J n ■ A(\n)d\ 



(12) 



This redefinition has the effect of decoupling usoft lines 
from collinear lines in the action, at the cost of introduc- 
ing Y factors into the operator Oq. However, we can 
see that since the Y are usoft and carry no large light- 
cone momenta they will cancel in Oq as a consequence of 



5 If we were to treat the photon like a hadron, then the analysis is 
almost identical to the case of the EM form factor discussed in 
later paragraphs. 



unitarity. The F's will however, show up in sub-leading 
Lagrangian and external operators. I will drop the hatted 
symbol from here on, and thus the reader should assume 
that all fields have been redefined as in Qllfl. 

The matching onto SCETII is trivial at this order. The 
off-shellness of the external collinear lines is reduced to 
being less than AQ. Taking the matrix element of this 
operator between the vacuum and one pion state yields 
the usual leading order result in terms of the pion wave 
function J5J. 



Power Corrections 

Power corrections arise from either matching onto 
higher order operators or from including corrections from 
the sub-leading Lagrangian into time ordered products 



uphngs be- 
tween usoft and collinear fields and is given by [l2llia . ll^| 



(TOP's). The order A action introduces couplings be- 
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with V 1 = n f2 n-D c /2 + + n^n-D/2, g.f. denotes 
gauge fixing terms, and 

igB; ± = [in ■ D c , iD^ c ] = n u {G n y^ . (14) 

The collinear gauge invariant field strength is 



We will also need the order A 2 Lagrangian 



n ■ D, 



— - — n ■ D u — - — Et] ^£ n . 
n-D c n-Dc C j2^ 

(16) 



As discussed above, it is important to understand that 
time ordered products corrections correspond to pertur- 
bations of the states in the effective theory. That is, 
in the effective theory the pion state is not the physical 
pion state, it contains only collinear modes in SCETII. 
The true pion is an eigenstate of the full Hamiltonian, so 
including perturbations into the time ordered product ac- 
counts for this difference in a s yst ematic fashion. There 
is a direct analogy with HQET |l5( which is perhaps illu- 
minating. In the limit where the quark mass is taken to 
infinity, the B meson is completely static. It is true that 
in the physical meson the quark has some kinetic energy, 
but these effects can be included when considering power 
corrections, via time ordered products with sub-leading 
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operators. In this way we can build up the full meson 
state order by order in the inverse quark mass. For a 
pedagogical discussion on the subject see (lif. 

Matching onto SCETI we may generate order A oper- 
ators by inserting B c ^~ into 0'°), 

0« =in v e x ^ n W) u J l5 (W^BP ± W) u)2 (WH n ) U3 , 

(17) 

But this operator does not interpolate for the pion 6 . 
Operators with insertions of can be absorbed into 
OW. We may also consider quark mass effects, which 
may arise from either the expansion of the fields via 

or from including mass effects in perturbative matching. 
Using a simple spurion[{| analysis and the fact that there 
is only one possible non-vanishing Dirac structure in the 
effective theory which violates chirality, it is simple to 
show that matching cannot generate any 0(A ) opera- 
tors 7 which are linear in the quark mass to all orders 
in perturbation theory. However, the introduction of a 
quark mass generates a new O(A) operator 

(19) 

The spontaneous breaking of chiral symmetry would lead 
these contributions to be numerically enhanced. 

The possible order A 2 operators which are bi-linear in 
the collinear quarks are 

=ie^ n WU(W^n-DW% 2 ^ 5 (W^ n U (20) 



0^=ie^ n W) wl (WI?lW% 2 (WI?lW%ai5(W^ n ) u 

(21) 

where the appearance of the momentum subscript implies 
the existence of non-trivial Wilson coefficients which ac- 

(2) 

count for possible insertions of the operator n- D c . O b , 
is a representative of a class of operators with two trans- 
verse covariant derivatives acting in all possible ways. 
The Wilson coefficients of these operators will be related 
by reparameterization invariance|l3j. Other operators 
involving collinear field strength operators can be ex- 
pressed in terms of linear combinations of these oper- 
ators. In addition, it is possible to generate operators 
with one collinear and one usoft quark, i.e. £ n W$ 75<7 US , 
but the contribution from this operator will be order A 3 , 
since an additional insertion of a sub-leading operator is 



6 We work in a frame where the transverse momentum of the pion 
vanishes. 

7 Though the operator itself is leading order it would come in sup- 
pressed by rriq/Q. 



necessary to get a non- vanishing matrix element between 
pion and vacuum. 

We now match onto the lower theory SCETII, where 
external virtualities are restricted to be less than AQ. In 
doing so, the external states will pick out a subset of the 
collinear modes, whose transverse label momentum is of 
order A. In addition, the usoft modes get relabeled to be 
soft modes. For Wilson lines this transformation is de- 
noted by Y -> S. { a ] will now scale as 0(A 2 /Q 2 ), since 
n ■ D picks out the smallest component of the collinear 
field which now scales in this way. The operators in the 

(2) 

class of O b will have vanishing Wilson coefficients when 
matching onto SCETII since we set the external trans- 
verse momenta to zero. Formally this occurs because 

T^Z 11 = 0. (22) 

However, there are a class of identical operators with 
usoft derivatives as well which will give A 2 /Q 2 corrections 
in SCETII. This scaling arises because usoft modes in 
SCETI match onto soft modes in SCETII, which scale 
as A/Q. In addition, we will generate new operators by 
considering the time ordered products in SCETI 

T 2 = J d 4 xT(0^cf\xf) , (23) 

where and c\ 2 ^ are any of the first and second order 
Lagrangian corrections respectively. Since the external 
states have vanishing transverse momenta as far as the 
label operators in SCETI are concerned, any operator 
with an odd number of transverse derivatives vanishes. 
TOPs with an even number of derivatives need not van- 
ish since there exists non- vanishing Wick contractions 8 . 
The TOPs T\ and T 2 however, will have non-vanishing 
matching coefficients and will contribute at order A/Q. 
To see how this scaling comes about, note that while 

(1 2) 

the power corrections to the Lagrangian C\ ' will scale 
down from being order A/Q to A 2 /Q 2 , this will not hap- 
pen within the time ordered product. The reason for this 
is that once we lower the virtuality of the external lines 
down to q 2 ~ A 2 , the collinear field scaling goes from 
\/AQ to A, but internal collinear lines connecting opera- 
tors at distinct points will still scale as they did in SCETI. 
This leads to enhancements in time ordered products rel- 
ative to "local" operators. For instance, if we consider 
T 2 , then a Wick contraction of collinear lines between the 
two operators leads to a fermion line which is off-shell by 
AQ. Thus the effective scaling of the product of these two 
fields is A/Q, whereas if they carried offshellness A 2 , the 



Note this is not true for the "local" operators OY 



scaling would be just A 2 /Q 2 , hence the enhancement. 
Furthermore, if the scale a/AQ is perturbative we may 
integrate out these off-shell modes perturbatively. 

Thus the correction to the leading order result may be 
written formally as 

Q 3 



'(1) 



^6F^=ld*xd*yl ^C{^ n , p \O^Lf\x)Cf\y)\Q) 



+ d A x j dujC{w^ n9 \O^Cf\xP). (24) 



Note that the large light cone momenta (u>i) flowing out 
of the leading order operator will now flow into the ver- 
tices of the sub-leading Lagrangian insertions. 

It is interesting to note that in the SCET formalism 
all the power corrections are factorable, in the sense that 
they can be written as products of matrix elements of var- 
ious types of fields 9 . That is, factorization is manifest 
in all the power corrections, since after the field redefi- 
nition, none of the various types of fields communicate. 
As a consequence of this, all of the soft fields may be 
factored into vacuum matrix elements which are inde- 
pendent of the hadron, since these fields have no overlap 
with the hadrons. This implies that, under our working 
assumption that SCET, in particular SCETII, as formu- 
lated in 0,^3 i s the appropriate effective field theory for 
the above processes, the soft pieces of all hadrons are uni- 
versal. Thus, there is hope that we can extract the soft 
structure functions and use them to make predictions in 
disparate processes. 



THE ELECTROMAGNETIC FORM FACTOR 

Leading Order 

Let us now consider the EM form factor. As we will 
see, the existence of two jets will have important ram- 
ifications. In this case, the Lagrangian splits into two 
pieces, one for each type of collincar mode, which do not 
communicate, 



(25) 



The usual leading order Brodsky-Lepage (BL) result 
was regained in SCET by matching the full QCD current 
onto four quark operators which are generated at order 
a.(Q)0, 

= 



(n + n) p x 




FIG. 1: Integrating out modes with virtuality yJAQ lead to 
operators which represent the "soft" piece of the wave func- 
tion. The solid (hatched) circle represents modes off-shell by 
Q W*Q)- 

The subscripts of these operators denote their color rep- 
resentation. 

Scaling the fields by the usoft Wilson lines has no ef- 
fect. The fields of opposing directions do not couple and 
therefore, the matrix element of these operators between 
pion states of opposing light like directions factorizes into 
the result (JSJ, with the octet contribution vanishing. 



Power Corrections 

When we're concerned with power corrections, we must 
consider the matching onto operators which may not nec- 
essarily have non- vanishing overlap with the pions. A 
non-zero overlap can be achieved via some TOPs' with 
subleading terms in the Lagrangian. In the EM case we 
may also match onto a quark bilinear 



(27) 



where the possible Dirac structures, T are given in JSJ. 
This operator is enhanced by a factor A -2 relative to the 
four quark operator since each collinear fields scales as A. 
However, charge conjugation implies that the only struc- 
ture with non- vanishing Wilson coefficient is 7^, which 
will not contribute for the case of the pion but will for 
the case of the p. At next order in the matching the 
operators 



o 



(-1) 



(28) 



(29) 



O 



(0) 



(n + n) p x 



are generated at order a®. Current conservation and RPI 
will relate the Wilson coefficients (C a ,Cb) of these two 
operators such that uJiC a + W4C& — 0. Note that we have 
introduced here two covariant derivatives with labels n 
and n, referring to the two distinct collinear sectors. 
(26) As opposed to the previous sub-leading operators, O a .b 
are not invariant under the ultra-soft field redefinitions 
In the redefined basis we have 



( - 1] =%(6ife WnY n ) Ul (Y^W^ n ) u 



(30) 



I thank Mark Wise and Iain Stewart for emphasizing this point . _ > 

tome. Of 1 ' = n^CnWnY^O^tW^PinCn)^ (31) 
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The leading order matrix elements between back to back 
pions of this operator vanish. To generate a non-zero 
overlap, all that is needed is the proper insertion of sub- 
leading operators which will inject one collinear quark 
into each jet. Furthermore, to get a non-vanishing Wil- 
son coefficient the TOP should contain an even number 
of insertions of This can be accomplished via 

a usoft (which become soft in SCETII) partonic fluctua- 
tion. The TOP 

T us (ui) = fd^XiTHCaicJu^Oi-^q^ixt) 

(32) 

gives an order A/Q contributions to the EM pion form 
factor of the form 



Q 2 SF 7T 



du>i(-K n , p | T us {uji) + h.c. | itn,p')- (33) 



This type of corrections was anticipated in [2j and is sim- 
ply a correction to the pion state. Note that each of 
the operators must be accompanied by a subleading La- 
grangian which contains a collinear gluon moving in the 
opposite direction relative the the gluon in the opera- 
tors. This insertion is needed to ensure that there are an 
even number of transverse covariant derivatives in both 
directions. 

If the scale AQ were perturbative then we could inte- 
grate the intermediate gluons as depicted in figure 1 . The 
ultra-soft quark lines, being fluctuations of virtuality A 2 , 
would only be closed once the matrix element is taken. 
There is an additional operator which gives an identical 
contribution but with n <-> n. TOPs at this order in- 
volving insertions of the quark mass have zero matching 
coefficients as they involve powers of the external trans- 
verse momentum. However, there will be quark mass 
dependence in the non-perturbative matrix elements in 
SCETII which will be relevant for 577(3) H3 and are only 
down by m/A. 

Notice that the usual leading order BL result is pro- 
portional to a s (Q) while this soft contribution scales as 
a s (y/QA) 2 . The data only reaches Q 2 = 10 GeV 2 , with 
error bars as large as the signal itself at larger Q[l^|. 
Given that A is of order 1 GeV, the scale y/AQ is likely 
non-perturbative over most of the range of the data, giv- 
ing the TOP an enhancement of l/a s (Q) relative to the 
leading order BL result. Thus, we see that the lack of 
concordance between theory and data may be due to 
the enhanced power correction in the EM form factor. 
Of course, there are other possible reasons for the dis- 
crepancy. It could be that using the asymptotic wave 



function is a poor approximation for these values of Q. 
However, the fact that the transition form factor seems 
to agree with the data0]> within theory errors, lends 
credence to the possibility that the enhanced power cor- 
rection discussed here could be the real culprit. Finally, 
the discrepancy could also be due to the extrapolation of 
the 7*p — * irn data to the pion polefl9|. 
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